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Abstract 

The analytic continuation of the Melhn transforms to complex values of N for the basic 
functions gi (x) of the momentum fraction x emerging in the quantities of massless QED and 
QCD up to two-loop order, as the unpolarized and polarized splitting functions, coefficient 
functions, and hard scattering cross sections for space- and time-like momentum transfer 
are evaluated. These Mellin transforms provide the analytic continuations of all finite 
harmonic sums up to the level of the threefold sums of transcendentality four, where the 
basis-set gi{x) consists of products of NiELSEN-integrals up to transcendentality four. The 
computer code ANCONT is provided. 
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Nature of the physical problem: 
Calculation of the analytic continuation of 
the Mellin transforms to complex values of N 
for the basic functions gi{x) emerging in the 
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QED and QCD up to two-loop order. 
Method of solution: 

Serial representations for complex argu- 
ments and approximate polynomial represen- 
tations. 

Restrictions on complexity of the problem: 
non 

Typical running time: Prom fractions of a 
second for individual representations to 3:30 
minutes for the test of all options on a HP- 
9000/755 workstation. 



1 



1 Introduction 



Hard processes in massless gauge field theories as QED and QCD exhibit MELLlN-structures. 
The observables a can be written as Mellin convolutions 

a{x) = \A®B\{x) = / dxi / dx2 5{x — xiX2)Ai{xi)A2{x2) ■ (1) 

Jo Jo 

Here, the functions Ai{xi) denote hard cross sections and splitting functions, respectively. In 
practice even multiple Mellin convolutions occur in the calculations of higher order Feynman 
diagrams. Integrals of this type can be diagonalized by the Mellin transform 

M \a{x)] (iV) = dx x^a{x) = M Mxi)] (iV) ■ M Walxa)] (iV) (2) 
Jo 

and a{x) can be found by a single inverse Mellin transform after evaluating the functions 
M [Aj(xj)] (N), which are simpler in general. For e N"*" the Mellin transforms of the functions 
Ai{x) can be represented in terms of linear combinations of finite harmonic sums [1-3] and their 
polynomials. 

Q (m ^ i^^Mh)r ^ [^^Mh)r [sign(fc^)]"- 

ni=l Til 112=1 nm=l 

Up to the level of two-loop order only harmonic sums up to rank 4 contribute. The set of 
functions Ai{x) is formed by polynomials of Nielsen integrals 



SnA^) = ) ' .X, , / - l0g"-'(^) l0g^(l - ZX) . (4) 



'"•^^^^ ~ {n-l)\p\ Jo z 

As was shown in Ref. the representation of the Mellin transforms in terms of finite harmonic 
sums may be reduced significantly by applying algebraic relations, cf. also 0, between these 
sums which are implied by index permutation and decomposition. Up to the level of two-loop 
order only 25 basic functions remain out of which the coefficient functions and splitting functions 
can be assembled as Mellin polynomials. 

Any of the multiple Mellin convolutions discussed above can thus be traced back to poly- 
nomials of the basic functions and a single numerical Mellin inversion which is performed as 
a complex contour integral. The Mellin transforms, which are firstly obtained at the positive 
integers, have to be analytically continued to the complex A^-plane. This analytic continuation 
is unique 0. For single harmonic sums the continuation is well-known 

^,(A^) = (-l)'=-^^^^^('^-^niV + l) + 4 (5) 

S^,[N) = j^^P'^'-'KN + 1) - , (6) 

with ip{z) the logarithmic derivative of the Euler F-function and 

1 



f3{z 



(7) 



2 

ct = IE (8) 

4 = Cik), k>2 (9) 

Cl = log(2) (10) 

4 = (l-^)ak), k>2. (11) 
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Here 7^ denotes the Euler-Mascheroni constant and ({k) is the Riemann (-function. 

It is the aim of the present paper to calculate the analytic continuations of the Mellin 
transforms of these 25 basic functions and to provide a code for numerical evaluations. The 
paper is organized as follows. The basic functions are introduced in section 2 and their structure 
is discussed. In section 3 the representation of the basic functions is given for positive integers. 
These representations which are given in terms of multiple harmonic sums do also provide first 
detailed numerical tests for the representations used for the analytic continuations later. The 
Mellin transforms for complex argument are given in section 4 for the functions of the type 
+ x) and for the functions {fk{x) — /fc(l))/(x — 1) in section 5. Here we aim on high 
numerical accuracy. Section 6 describes the basic options of the code ANCONT. The results are 
summarized in section 7. 

2 Basic Functions 

Using the algebraic relations given in Ref. one may show that the linear representations up to 
threefold harmonic sums up to transcendentality four can be expressed by the single harmonic 
sums S±k{N) and the MELLIN-transforms of the following 25 basic functions : 

log(l + x) 
X + 1 

I N log^(l + x) 

X + 1 

93{X) = —— ^ 5-2,1 AT 
X + 1 

9.{x) = - S,MN) 

X + 1 

log(a:)Li2(x) 

X + 1 
X +1 

X + 1 

X + 1 

9io{x) = ^ ^ 5_i,_2,-l(iV),52,-l,l(iV) 

I + X 

log(l -x)Li2(x) 

9ll{x) = ^ 5_i,2,l(A^), 5-2,1,1 

gu{x) = log(l-a:) ^''^~f^ ^ 5_i,_2,-i(iV), 52.-i,i(iV), 5'_2,-i,-i(iV) 

X + 1 

^13(^) = ^^^^Li2(-X) ^ 5i,2,-l(iV),52,l,_l(iV) 

log^(l + x)-log^(2) 
^14(3;) = — ^ 5_i,i,_i(A^) 



(12) 
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^^^^^^ ^ Ml+fI_M2)Li2(x) ^ 5-i,-2,i(iV),52,-i,-i(iV),5-2,-i,i(iV) 

X — 1 

log(l + x) - log(2) 
fi-iela;) = Li2(-x) ^ 5_i,2 -i(iV), 5_2,i -i(A^) 

X — 1 

( \ log(a:) log^(l + x) 

X — 1 

Li2(x)-C(2) 

fi'isia;) = ^ ^ i'2,i(iV) 

X — 1 

Li2(-x) + C(2)/2 
^i9(a;) = ^ ^_2,-i(A^) 

X — 1 
X — 1 

X — 1 

^ /■W + (5/8)C(3) ^ S._,_,(A').S_,_„(7V) 

X — 1 



with 



h{x) = / — log(l - z) log(l + z) . (13) 

Here we listed also all harmonic sums of rank 4 which contribute to the respective Mellin 
transforms. 

The functions gi{x) belong to the class of NiELSEN-integrals |^ Sn,p(x), Eq. (p. The degree 
of trans cendentality of these functions is defined to be r = p + n. The transcendentality of the 
product of two of these these functions is the sum of their transcendentalities. The measures 
dx/{x ± 1) are defined to be of transcendentality 1. The polylogarithms are given by 

Li„(x) = ^pM. ^ s„_i,i(x) = j^ll^ f-\og-\z) log(l - zx) for n > 2 . (14) 
a log(xj (n — 2j! Jo z 

The logarithms log(l ± x) are related to the dilogarithm by [0 

(iLi2(±x) 



(ilog(x) 

and 

Similarly, 



Lii(±x) = -log(l Tx) (15) 



Lio(x) = . (16) 

1 — X 



holds. As a generalization of the NiELSEN-integrals (^) one may wish to consider 

Sn,p,,(x) = ; . . / - log""^(z) log^(l - zx) \og%l + ZX) . (18) 

[n — ly.plql Jo z 

Both the classes of functions (^ and (|18D are non-minimal and large subsets of them can be 
represented already by the polylogarithms Lii{y), where y denotes a function of x. Relations 
between the different NiELSEN-integrals and their generalizations are, however, more easily 
established using non-minimal representations, since the argument structure turns out to be 
simple in the latter case. Following this line we have introduced in the above set of functions 

h{x) = 5i,i,i(x) = / - log(l - z) log(l + z) = -Si,2(x2) - Si,2(x) - Si,2(-x), (19) 
Jo z 2 

instead referring to Si^2{x'^)- 

The Mellin transforms of the functions (|12D are related to a series of finite harmonic sums 
through which all harmonic sums up to level three and transcendentality four may be expressed 
in terms of polynomials, except of the well-known case of single harmonic sums and polynomials 
which are made of only single harmonic sums. To establish the connection more closely we 
mentioned above those sums in which the respective Mellin transforms occur in explicit form, 
see also Ref. 0. 



3 Representations of the Mellin Transforms at Positive 
Integers 



To compare the numerical expressions of the analytic continuations of the MELLlN-transforms 
given in the forthcoming sections we summarize the exact representations of the Mellin- 
transforms for e N"*" for the basic functions in terms of harmonic sums. These relations 
are used as one test for the numerical accuracy of the analytic continuations. 



M[g,{x)] (iV) 



1 



5_i,i(iV) + - log^(2) - [5i(iV) - 5_i(iV)] log(2) 



(20) 



M[g2{x)] (N) 



M[gs{x)] (N) 
M[g,ix)] (AT) 



M[g,{x)] (N) 



2i-iy' 5i,i,_i(Ar) - log(2) [5i,_i(A^) - 5i,i(Ar)] 



-^log2(2) [5i(Ar)-5_i(Ar)] + ilog3(2) 



\N+1 



\N+1 



S_2,i(A^) - C(2)5_i(Ar) + -C(3) - C(2) log(2) 



S2MN) + log(2) [^2(A^) - S^2iN)] + -C(2)S_i(Ar) 



4c(3) + ^C(2)log(2) 



:-i 



S^2AN) + 25_3,i(A^) - 2C(2)5_2(A^) - ^C(2)' 



(21) 
(22) 



(23) 
(24) 
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M[ge{x)] {N) 



M[gr{x)] (iV) 



M[g,o{x)] (N) 



-1) 



N 



5_3,i(7V) - C(2)^-2(iV) + mS-i{N) 



+^C(2)^ - 2Li4 Q - ^C(3)log(2) + ^C(2)log^(2) - ^ V(2) 
(-1)^' js3,-i(iV) + log(2) [Ss{N) - S.,{N)] + ic(2)^-2(iV) 
-^C(3)5-i(iV) + ^C(2)^-|C(3)log(2)| 



(25) 



(26) 



M[gs{x)]{N) = (-1) 



N+l 



S-2,iAN) - C{^)S-i{N) + Lu Q - lc{2f - ^C(3) log(2) 



-ic(2)log^(2) + ^log^(2) 



M[g,{x)]{N) = (-l)^+^|52,i,-i(A^)+log(2)[52,i(A^)-52,-i(A^)] 
-^log2(2) [S,{N) - S.,{N)] - ^C(3)^-i(A^) 
-3Li4 Q + ^C(2)^ - ^C(3)log(2) + ^C(2)log^(2) - ^log^(2)} 
{-lf+'I^S_2,-i,-m + S2,-iAN) - log(2) [S.^AN) - S.2,-im 



+ 2 



C(2) - V(2) [S2{N) - S^2{N)] + -C(3) [5_i(7V) + log(2)] 



(27) 



(28) 



(29) 



M[gu{x)]{N) 



M[gu{x)]{N) 



(-l)^|5_i,2,i(7V) + 2S.2,i,i{N) - 2C(3)5_i(Ar) - C(2)5_i,i(iV) 
-gc(2)^ + 3Li4 Q - ic(2)log^(2) + iV(2)| 

(-l)^|52,-M(iV) + ^-2,-l,-l(iV) + ^-l,-2,-l(iV) 



(30) 



- log(2) [S-2,i{N) + S-i,2m + -C(2)5_i,i(7V) 

+ [S_,{N)S_2{N) + Ss{N)]\og{2) 



+ 2 



C(2) - log2(2) [S^iN] - S.^m + ^mS-,{N) 



-4Li4 ( -) + lo^r - |C(3) log(2) + ^C(2) log^(2) - ^ log^(2) I (31) 



M[g,s{x)] (TV) 



(-l)^|5i,2,-i(iV) + 2S2,rMN) + -C(2)5i,-i(iV) 
-log2(2)[52(iV)-^_2(iV)] 

+ log(2) \S2,i{N) - Si,_2{N) - 2S2,-i{N) + Si{N)S2{N) 
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M[g,,{x)]{N) 



M[g,,{x)] (iV) 



M[g,,{x)]{N) 



+S3{N) - -C(2)^-i(iV) 
-3Li4 Q + la2r - f C(3)log(2) + ^C(2) W(2) - ^log^(2)| (32) 



M [gu{x)] (AT) = 2S_,,^,_^{N) - 21og(2) [S_,,_,{N) - S_^,,{N)] - log'{2)S_,{N) 



ic(3) + C(2)log(2)-^V(2) 



(33) 



5_i,_2,i(iV) + 52,_i,_i(iV) + 5_2,-i,i(iV) - C(2)5_i,_i(iV) + -C(3)5_i(iV) 
^C(3) - C(2) log(2)] [S,{N) - 5_i(7V)] + log(2) [S^^N) - 2C(3)] 
(C(2) - log2(2)) [S2{N) - S.2{N)] - log(2)C(2)5i(7V) 
+XC(2)' 



Lid 



^C(3)log(2)-ic(2)W(2)-^log^(2) 



2^_2,l,_l(iV) + S_i,2,-l(iV) 

+ l0g(2) [2 (5_2,l(iV) - 5-2,-l(iV)) + ^_l,2(iV) - ^-l,-2(iV)] 

+ ^C(2)^-i,-i(iV) - log2(2) [S.,{N) - S^iN)] - lmS^{N) 
- {\m - \m log(2)) [S.,iN) - S^iN)] 
+ log(2) [s^2,-i{N) - log(2) {S2{N) - 5_2(iV)) + ^C(2)^i(A^) 
-|^C(2)^ + 4Li4 (-) + ^C(3) log(2) - |C(2) log^(2) + - log^(2) 

-2 [5_i,2,-l(iV) + 5_i,i,_2(iV) + 5_2,l,-l(iV)] 

+ l0g(2) [5-l,-2(iV) - ^_i,2(iV) - 5_2,l(iV) + 5-2,-l(iV)] 

+^C(3)^-i(A^) + ^W(2) [S_,{N) - S,{N)] - lc(2)5-i,i(A^) 
- 4Li4 (^) - ^C(3) log(2) + ^C(2) log^(2) - i log^(2) 



M[^i8(x) 




= -52,i(A^) + 2C(3) 




M[^i9(x) 


\iN) 


= -,5_2,-i(iV)+log(2)[52(A^)- 


'5-2(A^)] - ^C(3) 


M[g2oix) 


]iN) 


= 53,i(iV)-C(2)52(iV) + ic(2)^ 




M[g2i{x) 


]{N) 


= -52,i,i(A^) + ^C(2)^ 




M[g22{x) 


]iN) 


= 253,i(iV) + ^54(A^) + i52'(A^) 


- 2C(2)52(A^) + ^C(2)^ 


M[g23{x) 


]{N) 


= ,S_3,_i(iV)-log(2)[,53(7V)-,5 


-3(A^)] + k{2)S2iN) 



M[g2,{x)]{N) = -S2,-i,-i{N)-S^2,-i,i{N)+log{2)[S2,i{N)-S2,-i{N)] 



+- 



C(2) - V(2) [52(A^) - S^2m - -C(3)5i(iV) 



(34) 



(35) 



(36) 
(37) 
(38) 

(39) 

(40) 

(41) 
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+ic(2)^ - 2Li4 Q - Im log(2) + lc(2) log^(2) - ^ log^(2) (43) 
M[(725(x)](iV) = -5_2,i,-i(iV)-log(2)[5_2,i(iV)-5_2,-i(iV)] 

+ ^log^(2) [5_2(iV) - S,{N)] + ^a^r (44) 

The above relations can be checked by numerical integration directly. In the code ANCONT we 
use the routines DAIND, Ref. [Q, and verified these relations numerically for the integer moments 
e [1, 20] at an accuracy of better than at least 4 ■ 10~^, and in many cases even of 0(10^^^). 



4 Mellin Transforms for the Functions f(x)/(x + 1) 



For functions of the type 



gi{x) 



x + 1 



and fi{x) e C°°[0, 1[, cf. 0, one may expand fi{x) into a Taylor series around x = 0, 

oo 
k=0 

The MELLlN-transform is then given by 



M 



x+1 



{N) = Y,akf3{N + k + l) . 



(45) 



(46) 



(47) 



fc=0 



Since the MELLlN-transforms of the type M[log'"(x)/(a;)](A^) can be calculated using the relation 

QN 



M[log'"(x)/(x)] (iV) 



M[/(x)] (iV) 



one obtains 



M 



iog™(x; 



l + x 



k=0 



The function (z) obeys the recursion relation 



/3M(2 + 1) = (-l)'»m! 



For the polylogarithms Lii(±x) one obtains 

Li;(±x 



M 



log'"(x) 



a; + 1 



oo (A_'\\k 



k=l 



(48) 
(49) 

(50) 

(51) 



as an example. Although these representations are quite general and may be applied to other 
NiELSEN-integrals as well, the corresponding series may not converge fast enough. 

Alternatively, the analytic continuation of the Mellin transforms containing the factor 
{x + 1)~^ can be performed using the transformation 



M 



l+x 



(AT)^/ dx^-f{x) = log(2) ■ /(I) - / dxx^-Mog(l + x)[Ar/(x)+x/'(x)] 

JO X + 1 JO 



log(2) ■ /(I) - E 4'^ {^^M[/(x)] (N + k-l) 



k=l 



+ M[f{x)]{N + k)} 



(52) 
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if f{x) is differentiable in ]0, 1[. The function log(l + x) can be approximated by 



in(i + x) ^ 



(53) 



k=l 



with an accuracy better than 3 x 10 ^. The polynomial is determined using the approximation 

(1) 



given by the MINIMAX-method with an adaptive choice of arguments x. The coefficients 
are given by 



4" 
4" 
4" 

,(1) 



0.999999974532238E + 
0.333203435557262E + 
0.191451164719161E + 
0.792107412244877E - 1 
0.538406208663153E - 2 



a?) 



-0.499995525889840E + 
-0.248529457782640E + 
-0.137466222728331E +0 
-0.301109656912626E - 1 



Similar representations are found in the literature cf. [11-13] and have been used in Ref. [|l| 
before. Up to two-loop order also the representations of log'^(l + x) are needed with A; = 2, 3, 



In^fl +x) 



11 

E 

k=2 

13 



\n^{l + x) ~ E4^^^^ 

A;=3 



(54) 
(55) 



The coefficients are 



and 



(2) 

a\ 

(2) 
(2) 
(2) 



a 



10 



(3) 

(3) 
«5 

(3) 
% 

(3) 
"11 
(3) 



a 



13 



0.999999980543793E + 
0.916516447393493E + a 
0.745873737923571E + 
0.467104011423750E + a. 
0.936814286867420E - 1 



(2) 



(2) 
(2) 
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hi 



0.999999989322696E + 
0.174988008499745E + 1 
0.191539974617231E + 1 
0.162987195424434E + 1 
0.628710122994999E + a 
0.328953352932140E - 1 



(3) 
a\ 

"6 

"8 

(3) 
"10 
(3) 
12 



-0.999995797779624E + 
-0.831229921350708E + 
-0.634523908078600E + 
-0.261348046799178E + 
-0.156249375012462E - 1 



-0.149999722020708E + 1 
-0.187296689068405E + 1 
-0.185963744001295E + 1 
-0.117982353224299E + 1 
-0.211307487211713E + 



Eqs. (|5^j55|) hold at an accuracy better than 2 x 10^^ for x e [0, 1]. 

The MELLlN-transforms of the functions log(l + x)/(l + x) and log^(l + x)/(l + x) are 
obtained by 



M 



log^(l +x) 
1 + x 



(AT) 



1 



k + l 
1 



log^+^(2) - A^M log''+'(l + x) (A^-l) 



fc+i. 



k 



iog'=+^(2) - iv y: 



Li(fc+1) (k+l) 



1 N + l 



(56) 
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with Li{2) = 11 and Li{3) = 13. Similarly the MELLlN-transform of the function gnlx) is 
obtained by 



M 



\og{x) 

X — 1 



log^(l + x) 



11 



(N) = E 4 V'(iv + k + i). 

k=2 



(57) 



The remaining MELLlN-transforms are obtained from the representation Eq. (|5^ ) and are ex- 
pressed in terms of the known MELLIN-transforms for the functions f{x) and f'{x), cf. : 



M[Li2(x)] (iV) 
M[LUx)]iN) 



1 



1 



C(2) 



Si{N + l) 
N + l 



M [log(l - x)] {N - 1) 



SijN) 
N 



etc. The single harmonic sums S±k{N) are given in Eqs. (|^,|). 
One obtains the representations : 



M 



M 



Li2(a:) 
1 + x 

Li2(-a;) 
1 + x 



(AT) 
(AT) 



C(2)log(2)-E 



k=l 



N k 



N + k 



--C(2)log(2) + 5: 



,(1) 



\ AT + 2 (AT + 



[log(2) 



-p{N + k + l)] 



M 



log(x 
M 



^LisOr) 
1 + x 

LigOr) 
1 + x 



(AT) 
(N) 



E 



.=1 (^^ + ky 

C(3) log(2) 



C(2) + V''(fc + 1) -2 



Si{N + k) 
N + k 



-E 



(1) 



^ C(3). ' 



k=l 



N + k 



{N + ky 



C(2)- 



SijN + k) 
N + k 



M 



Li3(-a^) 
1 + x 



(AT) 



--C(3)log(2) + E 



^?c(3)+^ic(2) 



k=l 



N + kA' 



{N + kY2 



M 



Si,2(a^) 
1 + X 



(AT) 



(AT + A;) 
log(2)C(3) - E 



3[log(2)-/5(Ar + fc + l)] 



k=l 



N 



N + k 



C(3) 



A; 



(A^ + fc)2 2 



Sl{N + k) + S2{N + k) 



M 



Si,2(-a;) 



1 + x 



(AT) 



-C(3)log(2)-E 



A^ 



fe=i 



N + k 



.(2) 



C(3) _ ^ 
8 2f^^N + k + l 



If 

2,4^3 AT + fc 



(58) 
(59) 



(60) 



(61) 
(62) 



(63) 



(64) 



(65) 



(66) 



The function 81,2(3^^) does not occur singly but rather in the combination Si,2(a;^)/2— Sl,2(2;) — 
Sl,2(— a^) in the MELLIN-transforms which are related to the harmonic sums up to the level being 



10 



investigated in the present paper. It emerges via the integral 



rx rl 1 

h{x) = / - log(l - z) log(l + Z) = -Si,2(x2 
JO z 2 



81,2(3;) - Si,2(-a;) . 



(67) 



We firstly calculate the MELLlN-transform for this integral. The MELLlN-transform for 81,2(2;^) 
is then easily obtained by a linear combination using the relations given above. 



M 



h{x)' 
l + x 



11 



(iV) 



--C(3)log(2) + ^ 



mSi{N + k) 



fc=l 



N 



k=2 
5 



N + k 



N + k 



c(3)-E 



1=1 



^)SiiN + k + l)' 
N + k + l 



(6J 



N 




1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 



Eq. (39) 



-0.717E-03 
-0.117E-02 
-0.160E-02 
-0.207E-02 
-0.257E-02 
-0.308E-02 
-0.359E-02 
-0.410E-02 
-0.461E-02 
-0.511E-02 
-0.559E-02 
-0.607E-02 
-0.652E-02 
-0.697E-02 
-0.740E-02 
-0.782E-02 
-0.823E-02 
-0.862E-02 
-0.900E-02 
-0.936E-02 
-0.972E-02 



Eq. (70) 



0.969E-04 
0.205E-03 
0.292E-03 
0.378E-03 
0.452E-03 
0.505E-03 
0.531E-03 
0.529E-03 
0.501E-03 
0.449E-03 
0.376E-03 
0.286E-03 
0.181E-03 
0.637E-04 
-0.630E-04 
-0.197E-03 
-0.338E-03 
-0.483E-03 
-0.632E-03 
-0.784E-03 
-0.936E-03 



Eq. (60) 



0.212E-08 
0.442E-08 
0.761E-08 
0.117E-07 
0.169E-07 
0.232E-07 
0.307E-07 
0.397E-07 
0.501E-07 
0.622E-07 
0.758E-07 
0.909E-07 
0.107E-06 
0.125E-06 
0.144E-06 
0.164E-06 
0.184E-06 
0.205E-06 
0.226E-06 
0.247E-06 
0.268E-06 



Table 1: Comparison of the relative accuracy of different approximations for M [Li2(a;)/(1 + x)] (N). 

Since the harmonic sum S'_2,i ( A^) contributes already to the next-to-leading order anomalous 
dimensions, approximations for M [Li2(x)/(1 + x)] (N) have been used in the literature for some 
time [n, 1^]. Two approximations used before are 



or 



M 



M 



LiaOr) 
l + x 

Li2(x) 



l + x 



(AT) 



(iV) 



1.01 



0.846 



1.155 



1.074 0.55 

+ 



+ 2 + 3 A^ + 4 A^ + 5 + 6 



1.004 



0.846 1.342 

+ 



1.532 0.839 

+ 



A^ + 2 A^ + 3 A^ + 4 N + 5 N + 6 



(69) 



(70) 
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cf. IIJ] , which can be compared to our representation Eq. (|60D . In Table 1 the integer moments of 
these representations are given for e [0, 20]. Eq. ( pDj ) is more precise than Eqs. ( |5P| , [7D| ) by three 
to four orders of magnitude. Other quantities were approximated by similar representations as 
dB^JTg) in ^ recently. 

A MELLlN-transform of a function carrying the factor log(l — x)/{l + x) can be performed 
if the remainder factor obeys a polynomial representation, since the former factor may be trans- 
formed into the MELLlN-transform of log(l + x)/{l + x), cf. Ref. 0. Eqs. (44,45), by 



log(l 



M 

1+x 
Using Eq. 
"log(l 



x] 



(N) 



M 



log(l + x) 



X 



{N)-(3{N + l)[ij{N 



7E-log(2)]+/5'(Ar+l). 

(71) 



M 



) one obtains 
x)Li2 



-x] 



1 + x 



(AT) 



E 

k=l 



,(1) 



11 



log2(2)-^a| 



(2) 



N + k 



1=2 



N + k + l 



-(3'{N + k + l) 



+P{N + k + l) [Si{N + k)- log(2)] 



(72) 



The MELLlN-transform of the function log(l — x)Li2(a;)/(l + x) can be calculated similarly 
to the case which was discussed before. Li2(x) may be represented by 



Li2(x) ~ Po^'Va;) + 



Pi^\x) + \og{l-x)PP{x] 



where Po^i2{x) are the polynomials 



,(1)/ 



11 



PI 



(1), 



k=0 

49 



(1), 



x] 



XI 



, 11 7 2 1 s 

h C(2) H X x^ + -x^ 

36 ^ 6 12 9 



1 3 
-X^ 

3 



with 



"0 
-2 

-6 
-8 
-10 



-0.283822933724932E + 
-0.124975762907682E + 1 
-0.280403220046588E - 1 
0.532318519269331E + 
0.138194913357518E + 1 
0.506649587198046E + 



Ci 
^11 



0.999994319023731E + 
0.607076808008983E + 
-0.181869786537805E + 
-0.107281686995035E + 1 
-0.111100841298484E + 1 
-0.100672390783659E + 



(73) 

(74) 

(75) 
(76) 



The representation Eq. ([75| ) holds at an accuracy of better than 3 x 10 for x e [0, 1]. 

The Mellin transform of log(l — a:;)Li2(x)/(l + x) is thus expressed as a sum over Mellin 
transforms M[log(l-x)/(l+a;)]A^), Eq. (|7T]) and M[log^(l-x)/(l+x)](Ar). The latter expression 
is related to the harmonic sum 5'_i^i^i(A^), 0, Eq. (60). This sum can be expressed by the sum 
'S'i,i,-i(A^) and lower order harmonic sums through 



S. 



,iN) = 5i,i,_i(Ar) + Si(Ar)5i,_i(Ar) + S„2,i(iV) + '5-i,2(iV) 



+ 2 



S',iN)S-,iN) - S^,iN)S2iN) -5_3(iV). 



(77) 
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The harmonic sum S'i^i^_i(A^) is related to M log^(l + + x) (A^), which can be calculated 
using the representation (^) and lower order Mellin transforms. A simpler representation of 



M 



log (1 — x)/(l + x) (N) is given by 



M 



log (1 — x) 
1 + 0; 



(N) 



E 



,(1) 



^."'^ {N + k 



Sl{N + k) + S2{N + k)\ - [Slik) + S2{k) 
1 



+ -C(3)-C(2)log(2) + -log^(2) 



(7J 



referring to (pSf). Similarly 
logfl — x) 



M 



(N) 



1 + x 

holds. The above relations yield 
log(l - x)Li2{x) 



E4 

fc=i 



(1) 



^5i(Ar + K)- S,{k)\ + \ [\og\2) - C(2)] . (79) 



M 



l + x 



(N) 



11 

k=0 



.(1) 



log2(2) - C(2) 



1=1 



I 



N + k + l 



S^{N + k + I) - Siil) 



+ E PSkila^) - C(2) log(2) + \ \og%2) 



i2 



N + k + l 

s'm-s2{i) 



Sl{N + k + l) + S2{N + k + l) 



(80) 



Finally, the MELLlN-transform of log(l + x)Li2(— x) /(l + x) is given by 



M 



log(l + x)Li2(-x) 
l + x 



(iV) ^ -lc(2)log^(2) + i{|-^ 



11 



+ E' 

k=2 



.(2) 



N 



N + k 



1 _ log(2) - P{N + k + 1] 
2^^ ^ N + k 



5 Mellin Transforms for the Functions f(x)/(x — 1 



For this class of MELLlN-transforms the role of log(l + x) in the previous section is taken by 
log(l — x), provided the functions /(x) in the numerator do vanish at a sufficient degree as 
X ^ 1. Unlike log(l + x), the function log(l — x) has no simple polynomial representation in 
the range x e [0, 1[ at a comparable accuracy to Eq. (|53D. One may be tempted to express the 
numerator functions in a series of log(l — x) instead. This is indeed possible for a wide class of 
Nielsen integrals as Si,p(x) in a simple manner, cf. also [l7, IS], 



^^•^^"^^ p\^^k\p + k 



\ogP+\l-x), 



(82) 
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with Bk the Bernoulli numbers.^ The Taylor expansion of log'^(l — x) is given by 

J 



loe'^d - x) 



1=0 L 



n 
m 



X 

J. 



(83) 



Here 
by 



n 
m 



denote the STIRLING-numbers of the first kind [|l^. The MELLIN-transform is given 



M [log^l - x)] (AT) = {-lfk\^S,,...,m = ^ B{N,b + 1)1,^0 



where B{a, b) denote Euler's Beta function. The finite multiple harmonic sums "S*!^ ...^ i(^) are 

given in f^, Eq. (158), and obey a determinant-representation. The elements of the correspond- 
ing matrix are numbers and the single harmonic sums 



^ 1 r_iv-i 



with ci = and q = ({I) for / > 1. In this way one obtains the analytic continuations. With 
growing values of k the explicit expressions become rather lengthly. For numerical computations 
one may use the recursion relation, Eq. (164), in Ref. 0, 



Si,...,i{N) = -J2SiiN)S,,...,,{N) 



1=1 



for complex values of A^. Alternatively to these exact expressions one might also use the Taylor 
series Eq. (^) for which the MELLIN-transform is easily calculated. The SxiRLlNG-numbers of 
the first kind can be represented using a recursion relation, cf. Eqs. (166,167) in Ref. @. The 
Mellin transform of Si.p(a;) is thus given by 



M [Si,,(x)] (AT) = ^ 



p! ^ k\ p + k N 



17) 



p+k 



On the other hand the representation 



C(n 1 



is obtained by partial integration. This relation is more compact than the former. Eqs. ( ^^8]) 
imply an interesting relation between the multiple finite harmonic sums of a single index and 
the Bernoulli numbers and the C^function. 

Since not all Nielsen integrals which we are considering in the present paper can be expressed 
easily as a series in log(l —x) and the corresponding series are found to converge not fast enough 
in a series of cases we are going to apply a somewhat modified representation. To guarantee a 

^For other Nielsen integrals, as already for Li;(a;), I > 3, the serial expansion in log(l — x) leads no longer to 
optimal representations due to the emergence of nested multiple series, see for an example. 
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fast convergence also in the range of large values of a; < 1 we subtract from some of the numerator 
functions a polynomial in log(l— a;) and x of low degree. The remainder function is expanded into 
a polynomial by the MINIMAX-method. In both cases the Mellin transforms can be calculated 
analytically afterwards. In some cases further partial integrations have to be performed. Let us 
now discuss the individual cases in detail. 

Some of the MELLlN-transforms are of the type 

^log*^(l + a;) -log'=(2) 



M 



X — 1 



(AT) 



(89) 



where k = 1,2. The first factor in Eq. (^) can be represented by a polynomial. Here the 
denominator is divided out of a polynomial representation of the numerator. The corresponding 
approximations are found applying the MINIMAX-method, 



log'(l + x) -log'(2) 



X — 1 



m 

k=0 



X 



(90) 



with 1/(1) = 8 and L{2) = 9. The expansion coefficients are given by 





= 


693147166991375E + 


4" = 


-0.306850436868254E + 


"2 


= 


193078041088284E + 


4" = 


-0.139403892894644E + 




= 


105269615988049E + 


4" = 


-0.746801353858524E - 1 




= 


427339135378207E - 1 


4" = 


-0.161809049989783E - 1 




= 


288664611077007E - 2 







and 



^(2) 

.(2) 
O4 

.(2) 
^(2) 



0.480453024731510E + 
-0.519463586324817E + 
-0.427765744446172E + 
-0.263827078164263E + 
-0.525105367350968E - 1 



J2) 
"3 



0.480450679641120E + 
0.479285947990175E + 
0.360855321373065E + 
0.146927719341510E + 
0.874144396622167E - 2 



The related MELLlN-transforms have the representation 



M 



log(l + x) - log(2) 



M 



X — 1 

log'(l + a;) -log2(2) 



M 



log(l + x) 



x — 1 
- log(2) 



X 



M 



log(l + x) 



1 

log(2) 



Li2(xl 



X 



Li2( 



-X] 



{N) 
(N) 
(N) 
(AT) 



8 

E 

fc=0 
9 

E 

fc=0 
8 

E 

fc=0 
8 

E 

fc=0 



N + k + 1 



l{2) 



N + k + 1 



N + k + 1 



C(2) 



SiiN + k + 1) 



N + k + 1 



:C(2) + 



N + k + 1 
log(2) -p{N + k + 2) 



(91) 
(92) 
(93) 



N + k + 1 



In the case of the remaining MELLlN-transforms of the type 

lXx)-f{iy 



M 



X — 1 



(N) 



(94) 



(95) 



15 



firstly the integral 

is evaluated leading to 

M 



F{x) 
/(^) - /(I) 



dz 



fiz) - /(I) 



X 



(N) = - NM [F(x)] (N-1) . 



(96) 



(97) 



We now seek for an appropriate representation of the functions F{x). The respective functions 
are : 



/ dz 
Jo 



L 



Li2(^)-C(2) 
z-l 

^ ^ u,{-z)+ cm/2 

z - 1 
Si,2(^)-C(3) 



[Li2(x)-C(2)]log(l-x) + 2Si,2(x) 



^^2{-x) + -C(2) 



I 

fX 

/ dz 
Jo 



10 z — 1 

Here one may use the representation 



log(l - x) + Ii{x) 
= [Li3(a:)-C(3)]log(l-a:) + ^Li^(a:) 
= [Si,2(a;)-C(3)]log(l-x)+3Si,3(a;) 



9 .(1) 



(98) 
(99) 
(100) 
(101) 

(102) 



We derive approximations for Li3(a;), 81,2(2;), hil{x) and Ii{x) to be able to perform the Mellin- 
transform using these semi-analytic expressions. Here it is of particular importance to account 
for the log'^(l — x) behavior as x — > 1. One obtains 

(103) 
(104) 
(105) 

where -Po,?,2(^) ^'^'S the polynomials 



Lis 


[x) 




[x) 


+ 


p(2) 


{x) + log(l 


— x] 




[x) 






Sl,2 


[x) 




{x) 


+ 


p(3) 


[x) + log(l 


— x] 


pP 


{x) 


+ log2 


[1-x)pP{x) 


Lil 


[x) 


^ Pt^ 


{x) 


+ 


p(4) 


[x) + log(l 


— x] 




{x) 


+ log2 


:i-x)pt\x) 



12 



(106) 



fe=0 



Pi%) ^ c(3)->) + i+(3C(2)-H),_(|c(2)-|).^.(ic(2)-^ 



X 



pP{x) 



-1 + -X - 2x' + -x^ , 



(107) 
(108) 



with 



^0 


= 0.480322239287459E + 


"-i — 


-0.168480825099837E + 1 


•-2 


= 0.209270571633447E + 1 


^(2) _ 


-0.101728150522737E + 1 


„(2) 
C4 


= 0.160180000661971E +0 


„(2) _ 

— 


-0.351982379713689E + 


^6 


= 0.141033369447519E + 1 


c^') - 


-0.353344124579927E + 1 




= 0.593934899678262E + 1 


c^') - 
— 


-0.660019998525006E + 1 


''lO 


= 0.466330491799074E + 1 


^(2) _ 


-0.189825521858848E + 1 


^12 


= 0.339773000152805E + 







16 



(3), 



X] 



Pf (x) 



k=0 

C(3) 
205 



2035 205 

1728 ^ 144^ 
25 23 



95 2 43 o 1 4 

X H X X 

288 432 64 



o 13 o 1 4 

X H X X H X 

144 12 24 36 16 



25 , 3 2 

h 2x x^ + 

24 2 



2 3 1 
-x^ ; 



(109) 

(110) 
(111) 
(112) 



with 



C4 

^6 



-0.243948949064443E - 1 


Ci — 


0.000005136294145E + 


0.249849075518710E + 




-0.498290708990997E + 


0.354866791547134E + 




-0.522116678353452E - 1 


-0.648354706049337E - 1 


J3) _ 
C7 — 


0.644165053822532E - 1 


-0.394927322542075E - 1 




0.100879370657869E - 1 



P}'\x) 



12 

E 

k=0 

257 205 



C(2) + C(2)^ 



144 72 

C(2) + 

36 18^^ " 



/167 25^,^,\ 



144 8 



, ^235 

V 36 



6<(2V + V 18 



K(2) a; 



35 
12 



26 



19 



14 



11 



— X H X X H X 

3 2 3 12 



V 24 12^^ ' ' 



40 A 2 /109 



;ii3) 



;ii4) 



3C(2) ) x' 



(115) 
(116) 



with 





0.192962504274437E + 


Ci 


= 0.000005641557253E + 




-0.196891075399448E + 1 




= 0.392919138747074E + 1 


C4 — 


-0.290306105685546E + 1 




= 0.992890266001707E + 


c^^) - 

^6 — 


-0.130026190226546E + 1 




= 0.341870577921103E + 1 


c^') - 

^8 — 


-0.576763902370864E + 1 




= 0.645554138192407E + 1 


J4) _ 
^10 ~ 


0.459405622046138E + 1 


'^11 


= 0.188510809558304E + 1 


^12 — 


-0.340476080290674E + 







These representations hold at an accuracy of better than 2 x 10"^ 3 x 10"^ and 2 x 10-^ 
respectively, for x e [0, 1]. 

In the case of the function /i(x), Eq. (p^Bf), the function log(l + x) in the integrand can be 
approximated by Eq. (|53|) . The integral can then be evaluated analytically and takes the form 



h{x) 



+ pP{x) log(l -x) 



;ii7) 
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with 

p!^'\^) = jlcfx' (118) 
fc=i 

P^\x) = jldfxK (119) 

fc=0 

The coefficients c^^'* and are algebraically related to the coefficients a^^\ Eq. (|53|). Their 
numerical values are 



Cl — 


-0.822467033400776E + 


^2 


= 0.887664705657325E 


- 1 


c^') - 

^3 — 


-0.241549406045162E 


- 1 


C4 


= 0.965074750946139E 


- 2 




-0.470587487919749E 


- 2 




= 0.246014308378549E 


- 2 


— 


-0.116431121874067E 


- 2 




= 0.395705193848026E 


- 3 




-0.664699010014505E 


-4 









"0 ~ 


-0.822467033400776E + 




= 0.999999974532241E + 


4' = 


-0.249997762945014E + 


df 


= 0.111067811851394E + 


rff = 


-0.621323644338330E - 1 


"5 


= 0.382902328987004E - 1 


<if = 


-0.229110370338977E - 1 


4" 


= 0.113158200819689E - 1 




-0.376387065979726E - 2 


4' 


= 0.598229109013054E - 3 



This representation holds at an accuracy of better than 3 x 10 * for x e [0, 1]. 
The MELLIN-transforms of the above polynomials are : 



.k=0 



M 



M 

m 



.k=0 



(iV) 
(N) 



m+l 



m A 

(-iy^'E]^^i,....i(iv + ^) 



fc=i 



(120) 
(121) 



Here, 



2! 
1 

3! 



Sl{N) + S2{N) 



Sl{N) + ?,Si{N)S2{N) + 2^3 (A^) 



;i22) 

:i23) 



The Mellin transforms of the functions gis{x) — (722(2;) are then given represented by 



M 



Li2(x) - C(2) 



X 



11 



(iV) 



Sl{N) + S2{N) -C(2)5i(iV) + 5: 



-.(1) 



fc=0 



k=0 



N + k 



Sl{N + k) + S2{N + k) 



M 



Li2(-x) + C(2)/2 



X 



(N) 



la2)S^iN)-j:^S^iN + k) 

k=l 



(124) 
(125) 



18 



M 



Usjx) - C(3) 
X — 1 



12 



(iV) 



.(2) 



Si{N + k) 



fc=0 



^i2(Ar + A;) + ^2(A^ + A;) 

2 



1 12 

1 — .(4) 



E 



AT 



M 



Si,2(x) - C(3) 

X — 1 



(iV) 



(4) gl(iV + fc) _ (4) ^KA^ + fc) + g2(iV + fc) 

2'^' N + k ^■'^ iV + A; 

1 



(126) 



-C(3)5i(iV) + ^ [5i'(iV) + 2S,{N)S2{N) + 2Ss{N) 



N 



[^i3(Ar + A;) + 3SiiN + k)S2iN + k) + 2SsiN + k) 



-Pr 



(3) 



2,fc 



Sl{N + k) + S2{N + k)]] 



(127) 



M 



log(x)Li2(x) 
a; — 1 



11 



(iV) 



^ci'V'(iV + fc + l) 



fc=i 



5^ P!^]l Sr{N + A;)7/;'(Ar + A; + 1) - -ij"{N + k + l) 
fc=o ' ^ 



(128) 



A set of other integrals we are going to re-write in terms of integrals of the type 
Jq dzf2{z)/{l + z), the MELLlN-transforms of which are then evaluated using Eq. (|52D. The 
following integral relations are obtained : 



dz 



Lisi-z) + (3/4)C(3) 
z-1 



U,i-x) + -C(3) 



-Li3(x) log(l + x) - dz— 
Jo 1 



log(l - x) + Li2(x) Li2(-x) 

Li3(^) 



+ z 



-^^ /i(^) + (5/8)C(3) 
z-1 



-^ /i.2(-^)-C(3)A 
z — 1 



h{x) + -C(3) 



log(l - x) - 2Si,2(x) log(l + x) 



+2 r dz^ 
Jo 1 



+ z 

1 



+- 



Si,2(-x) - -C(3) 
1 



log(l -x) - -Ii{x) log(l + x) 



2 Jo 1 + z 
Here one may represent Li3(— x) and Si^2(~a^) by 



Lu(-x) 



Si 2(-a;) 



9 a(^) 



fc=i 

^ fc=2 



(129) 

(130) 

(131) 

(132) 
(133) 



The Mellin transforms pf g22{x) — g25{x) are represented by 
^Li3(-x) + 3C(3)/3" 



M 



X — 1 



(N) = -lc(2)^ + C(3)log(2) + ^C(3)Si(Ar)-M 



Li3(x) 
1 + x 



(N) 
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(1) 



k=l 

C(2) 



N + k ^ 
+ Si{N + k) 



C(3)- 



C(2) 
N + k 

1 



1 1 

{N + ky k{N + k) 



(134) 



M 



/i(x)+5C(3)/8' 



X 



(iV) = -2C(3)log(2) + 2M 



(N) 



N 



fc=i 



2C(3) 



1 + x 

gi^(Ar + A;) + g2(A^ + fc) 
+ A; 



+ T.cf'j^^S,{N + k)-±d 



k=l 



N 



k=0 



N + k 



Sl{N + k) + S2{N + k) 



(135) 



M 



Si,2(-a^) - C(3)/8 
x-1 



(iV) 



^C(3)log(2) + lM 



/i(x) 
1 + x 



(iV) - -C(3)5i(iV) 



iV 



9 



Z=l 



(i)5i(iV + A; + 



;C(3) 
(136) 



N + k + l 

We finally would like to add a remark on a recent analysis of the two-loop coefficient functions 



of deep inelastic scattering pO|. The representation given there also contains the harmonic sum 



>S'-i,i,i,i(A^) which was not needed to express the individual Mellin transforms in Refs. ^1 



cf. Ref. 1^, which may be caused due to the algebraic relations being applied in f^. This sum 
has the representation 



S-i,i,iAN) = (-1) 



N- 



6^ 



M 



log^(l - x) 
1 + x 



(iV) - Li4 



(137) 



The Mellin transform of the function log^(l— x) / (1+x) for complex argument is easily obtained, 



M 



log^(l -x) 
1 + x 



(iV) 



fc=i 



.(1) 



\A^ + fc 



S^iN + k) + SSiiN + k)S2iN + k) + 2SsiN + k) 



Sl{k) + ?>Si{k)S2{k) + 2S^{k) 



6Li4 



;i38) 



The above relations allow to express the Mellin transforms of all functions emerging in the 
coefficient functions and anomalous dimensions of massless gauge theories up to two-loop order 
and to evaluate their analytic continuation from integer to complex arguments. 



6 The Code ANCONT 



6.1 General Structure 



The code ANCONT calculates the Mellin transforms of the basic functions gi{x) both for integer 
and complex values of the Mellin index N. The parameters of the code are initialized in 
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SUBROUTINE ACINI. 



The calculations are performed in 

SUBROUTINE ACRUN. 

The code ends with 

SUBROUTINE ACEND. 

6.2 USER routines 

The user may access parts of the code via the user routines UINIT, URUN and UOUT. These 
routines arc called in ACINI , ACRUN and ACEND, respectively. SUBROUTINE UINIT may be used to 
actualize the running parameters. Via URUN the user may built her/his own Mellin convolutions 
or other functions of the basic functions and their Mellin transforms for positive integer or 
complex Mellin index using the functions ACGi, FCTi or XCGi, which are described below, as 
a library. SUBROUTINE UOUT may be used as an output-interface at the end of the code. To 
transfer data between different user routines the user may define COMMON- blocks named 

COMMON/USxxxx/ . . . 

6.3 Initialization 

The main parameters and constants of the code are defined in SUBROUTINE ACINI. This routine 
calls the subroutines INVINI, DEFAUL and the user routine UINIT. SUBROUTINE INVINI sets 
parameters and constants related to the Mellin inversion. The default running parameters of 
the code are set in SUBROUTINE DEFAUL. These are : 

IRUN = 
ITESTl = 1 
ITEST2 = 1 
I TESTS = 1 
IMIN = 1 
IMAX = 26 
lAPP = 1 
NMIN = 1 
NMAX = 20 
EPS = l.OD-9 

The running parameters are printed by SUBROUTINE WROUT. The parameters ITESTi initialize 
tests of the code, which are inactive for ITESTi.NE.l. For ITESTl. EQ.l the representation 
of the moments of the basic functions gi{x) for positive integer index by harmonic sums are 
compared to those obtained by numerical integration. ITEST2 . EQ . 1 induces a comparison of the 
representations for the analytic continuation of the Mellin transform of the basic functions with 
those by the harmonic sums for positive integer index. For ITEST3 . EQ . 1 a comparison is made 
for the Mellin inversion using the analytic continuations of the Mellin transforms of the basic 
functions and the numerical representations of the basic functions in the range x e [10~^, 0.99]. 

IMIN and IMAX mark the index range of basic functions to be used. Likewise NMIN and NMAX 
set the minimum and maximum positive integer moment for the tests. 
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lAPP selects the representation of the analytic continuation for gslx) = Li2(x)/(1 + x) to 
compare the representations Eq. ( pDD lAPP = 1 and ( pD| ) IAPP=2, ( [7D| ) IAPP=3, respectively. 

EPS (REAL*8) is a pilot parameter for the numerical integration of the program DAIND p[ 
and denotes the relative numerical accuracy to be obtained. 

The above parameters are available through the COMMON-blocks 

COMMON/TEST/ ITESTl, ITEST2 , ITEST3 
COMMON/RUN / IRUN 
COMMON/MOMPA/ NMIN,NMAX 
COMMON/FUNPA/ IMIN,IMAX 
COMMON/IAPP/ lAPP 
COMMON/EP/ EPS. 



6.4 Running 

The code provides three main lines, which are initialized setting the parameters ITESTl, 
ITEST2, ITEST3 = 1, respectively. 

For ITEST1= 1 the representations of the integer moments of the basic functions gi{x) in terms 
of harmonic sums are tested comparing them with the corresponding numerical integrals. The 
value of the Mellin moment is calculated. The relative accuracy comparing both calculations 
RAT = VAL1/VAL2 -1 and the value of the moment are provided. The following test-output is 
obtained for k = 9, ggi^x) = Si,2(— + x), rii = 2,n2 = 20, val = l.OD-9. Here N denotes 
the Mellin index, RATk the relative accuracy comparing the numerical integration and the 
representation by harmonic sums, and VAL the value of the Mellin moment. 



N,RAT9,VAL= 


2 


-1, 


. 98330241 1 190380E- 


■12 


1, 


.784970126521851E- 


-02 


N,RAT9,VAL= 


3 


-1, 


.992739306899693E- 


■12 


1, 


.447619078834311E- 


-02 


N,RAT9,VAL= 


4 


-1, 


.742939126359033E- 


■12 


1, 


.216425178651454E- 


-02 


N,RAT9,VAL= 


5 


-1, 


.722177955798543E- 


■12 


1, 


.048399327332255E- 


-02 


N,RAT9,VAL= 


6 


-1, 


.469158128486470E- 


■12 


9, 


.208975593568651E- 


-03 


N,RAT9,VAL= 


7 


-1, 


.485589429250922E- 


■12 


8, 


.209008584718454E- 


-03 


N,RAT9,VAL= 


8 


-1, 


. 187272502534142E- 


■12 


7, 


.404084086538535E- 


-03 


N,RAT9,VAL= 


9 


-1, 


.256883486178140E- 


■12 


6, 


.742389451876160E- 


-03 


N,RAT9,VAL= 


10 


-1, 


.036948304999896E- 


■12 


6, 


. 188923213501677E- 


-03 


N,RAT9,VAL= 


11 


-8, 


.926193117986259E- 


■13 


5, 


.719203265306224E- 


-03 


N,RAT9,VAL= 


12 


-8, 


.733014311701481E- 


■13 


5, 


.315598831318883E- 


-03 


N,RAT9,VAL= 


13 


-7, 


.922551503725117E- 


■13 


4, 


.965094655720756E- 


-03 


N,RAT9,VAL= 


14 


-5, 


.925260282424460E- 


■13 


4, 


.657875778975199E- 


-03 


N,RAT9,VAL= 


15 


-6, 


.212808045802376E- 


■13 


4, 


.386402245697693E- 


-03 


N,RAT9,VAL= 


16 


-3, 


.986810881428937E- 


■13 


4, 


. 144786898528877E- 


-03 


N,RAT9,VAL= 


17 


-4, 


.581890422628021E- 


■13 


3, 


.928366556207149E- 


-03 


N,RAT9,VAL= 


18 


-1, 


.886268918838141E- 


■13 


3, 


.733399977739974E- 


-03 


N,RAT9,VAL= 


19 


-3, 


.899103262483550E- 


■13 


3, 


.556850973662370E- 


-03 


N,RAT9,VAL= 


20 


-1, 


.719735465144367E- 


■13 


3, 


.396229940556720E- 


-03 



For ITEST2= 1 the representations of the Mellin transforms of the basic functions gi{x) valid 
for complex arguments are compared to the representation in terms of harmonic sums at posi- 
tive integer argument. The value of the Mellin moment is calculated. The relative accuracy 
comparing both calculations RAT = VAL1/VAL2 -1 and the value of the moment are provided. 
The relative accuracy of these representations are given in Table 2 for = 2 to 20. 
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Table 2: Relative accuracy of the Mellin transforms comparing the approximative representations 
valid for complex arguments and the representation in terms of harmonic sums for e [2, 20] 
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882E- 


08 


-0 


248E- 


07 


8 


-0 


499E- 


08 





126E- 


07 





506E- 


08 


-0 


341E- 


09 





223E- 


07 


-0 


893E- 


08 





985E- 


08 


-0 


303E- 


07 


9 


-0 


525E- 


08 





157E- 


07 


-0 


377E- 


09 





612E- 


08 





257E- 


07 


-0 


963E- 


08 





109E- 


07 


-0 


365E- 


07 


10 


-0 


546E- 


08 





192E- 


07 





728E- 


08 


-0 


413E- 


09 





293E- 


07 


-0 


103E- 


07 





118E- 


07 


-0 


432E- 


07 


11 


-0 


562E- 


08 





231E^ 


07 





855E^ 


08 


-0 


449E^ 


09 





329E- 


07 


-0 


llOE- 


07 





128E- 


07 


-0 


506E- 


07 


12 


-0 


575E- 


08 





274E^ 


07 





993E^ 


08 


-0 


486E^ 


09 





366E- 


07 


-0 


116E- 


07 





138E- 


07 


-0 


587E- 


07 


13 


-0 


584E- 


08 





321E^ 


07 





114E^ 


07 


-0 


523E^ 


09 





404E- 


07 


-0 


123E- 


07 





147E- 


07 


-0 


675E- 


07 


14 


-0 


590E- 


08 





371E- 


07 





131E- 


07 


-0 


561E- 


09 





442E- 


07 


-0 


129E- 


07 





157E- 


07 


-0 


770E- 


07 


15 


-0 


592E- 


08 





426E- 


07 





148E- 


07 


-0 


600E- 


09 





481E- 


07 


-0 


135E- 


07 





166E- 


07 


-0 


872E- 


07 


16 


-0 


592E- 


08 





485E- 


07 





166E- 


07 


-0 


641E- 


09 





520E- 


07 


-0 


141E- 


07 





175E- 


07 


-0 


981E- 


07 


17 


-0 


589E- 


08 





549E- 


07 





186E- 


07 


-0 


683E- 


09 





560E- 


07 


-0 


147E- 


07 





184E- 


07 


-0 


llOE- 


06 


18 


-0 


583E- 


08 





617E^ 


07 





207E^ 


07 


-0 


729E^ 


09 





601E- 


07 


-0 


153E- 


07 





192E- 


07 


-0 


122E- 


06 


19 


-0 


576E- 


08 





690E^ 


07 





229E^ 


07 


-0 


776E^ 


09 





642E- 


07 


-0 


159E- 


07 





201E- 


07 


-0 


135E- 


06 


20 


-0 


.',()()E- 


0<S 





7(i9E- 


07 





2.11E- 


07 


-0 


S27E- 


09 





(i.S.'-SE- 


07 


-0 


1().')E- 


07 





21()E- 


07 


-0 


149E- 


0(i 



Table 2 : continued 
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For ITEST3= 1 the Mellin inversion for the basic functions gi{x) based on the representations 
of the Mellin transforms for complex argument are compared to the numerical expressions for 
the functions gi{x) in the range x e [10"'', 0.99]. The value of Qi as a function of x is calculated. 
The relative accuracy comparing both calculations RAT = VAL1/VAL2 -1 is given. The inverse 
Mellin transform to x-space is performed by a contour integral numerically. The singularities of 
the Mellin transforms for N e C are situated at the non-positive integers for the coefficient and 
splitting functions in massless fixed-order perturbation theory. The non-perturbative input 
densities are usually expressed in terms of polynomials of the type 

h{x) = J2Aix''^{l-xf^ , (139) 

i 

the Mellin transform of which is a sum of EuLER Beta-functions. Their singularities are as 
well situated on the real axis left of an upper bound. The inverse Mellin transform is given by 

hix) = / dz\m e**x-"(")/[c(^)] (140) 

JO L J 

where c{z) = co + ze**, cf. also |]l4i. The parameter integral over z can be performed by standard 



algorithms. The code segments the integration path logarithmically into N=20 or N=50 pieces on 
each of which the integral is calculated by the 8-point or 32-point Gauss formula. The starting 
point C is varied according to the rightmost singularity of the function to be inverted and is 
chosen close to it on the right. The angle of the linear path w.r.t. the real axis is chosen by 
$ = (3/4)7r for x < 0.8, $ = (7/8)7r for x > 0.8. and $ = (19/20)7r for x > 0.98. 

As an example the relative accuracy of the representation for the function giQ = 
[log(l + x) — log(2)] /{x — l)Li2(— x) between x = 10"'^ and x = 0.99 is 



X,RAT16,VAL= 


.lOOOOD-06 


-.44220D- 


08 




69315D-07 


X,RAT16,VAL= 


. lOOOOD-05 


-.10120D- 


07 




69315D-06 


X,RAT16,VAL= 


.lOOOOD-04 


-.13908D- 


07 




69314D-05 


X,RAT16,VAL= 


. lOOOOD-03 


-. 16350D- 


07 




69310D-04 


X,RAT16,VAL= 


.lOOOOD-02 


-. 15094D- 


07 




69267D-03 


X,RAT16,VAL= 


.lOOOOD-01 


.61713D- 


08 




68838D-02 


X,RAT16,VAL= 


.50000D-01 


. 94884D- 


08 




33499D-01 


X,RAT16,VAL= 


.lOOOOD+00 


-.19742D- 


07 




64836D-01 


X,RAT16,VAL= 


. 15000D+00 


-.92017D- 


08 




94219D-01 


X,RAT16,VAL= 


.20000D+00 


. 14400D- 


07 




12183D+00 


X,RAT16,VAL= 


.25000D+00 


.20299D- 


07 




14783D+00 


X,RAT16,VAL= 


.30000D+00 


.50955D- 


08 




17236D+00 


X,RAT16,VAL= 


.35000D+00 


-.15289D- 


07 




19554D+00 


X,RAT16,VAL= 


.40000D+00 


-.23563D- 


07 




21747D+00 


X,RAT16,VAL= 


.45000D+00 


-.13795D- 


07 




23827D+00 


X,RAT16,VAL= 


. 50000D+00 


.61787D- 


08 




25800D+00 


X,RAT16,VAL= 


.55000D+00 


.21289D- 


07 




27676D+00 


X,RAT16,VAL= 


.60000D+00 


. 19772D- 


07 




29461D+00 


X,RAT16,VAL= 


.65000D+00 


. 16678D- 


08 




31161D+00 


X,RAT16,VAL= 


.70000D+00 


-.19671D- 


07 




32783D+00 


X,RAT16,VAL= 


.75000D+00 


.22255D- 


05 




34332D+00 


X,RAT16,VAL= 


.80000D+00 


.22936D- 


05 




35811D+00 


X,RAT16,VAL= 


.85000D+00 


.30694D- 


05 




37226D+00 



^Some all order resummations, as e.g. the small-a; resumniation, lead also to singularities outside the real 
axis, cf. ||. 
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X,RAT16,VAL= . 90000D+00 
X,RAT16,VAL= .95000D+00 
X,RAT16,VAL= . 99000D+00 



.38631D-05 
-.25871D-02 
-.20484D-01 



-.38581D+00 
-.39879D+00 
-.40879D+00 



Similar accuracies are obtained also in the other cases. We would like to mention that for 
i = 3 and IAPP=2,3, respectively the accuracy is lower by three orders of magnitude up to 
X ~ 0.6 relative to the results for IAPP=1 and remains less above. 

The central functions of the code are : FCTi(N), XCGi, ACGi and FKNi The REAL* 8 FUNCTIONS 
FCTiCN), i=l . .26 provide the positive integer moments of the basic functions gi{x), N > 1 
obtained by numerical integration using the program DAIND [^. The integrands are defined in 
the REAL*8 FUNCTIONS FKTi(N) , i=l. .26. 

A second representation of the positive integer moments of the basic functions is given using 
the finite harmonic sums in the REAL*8 FUNCTIONS XCGi(N) , i=l. .26. 

The analytic continuations of the Mellin transforms are given in the COMPLEX* 16 
FUNCTIONS ACGi (N) , i=l . . 26. 

The basic functions gi{x) are given by REAL*8 FUNCTION FKNi(X) , i=l. .26. 

One may use the above options without performing the aforementioned comparisons. The 
corresponding pilot parameter is IRUN, which takes values between and 5. 

IRUN = : implies a test run. 

IRUN = 1 : The Mellin moments of the basic functions are calculated by numerical integration 
for positive Mellin indices. 

IRUN = 2 : The Mellin moments of the basic functions are calculated by their representation 
in terms of harmonic sums for positive Mellin indices. 

IRUN = 3 : The Mellin moments of the basic functions are calculated by their representation 
in terms of the analytic continuations ACGi for complex argument at positive integer moments. 

IRUN = 4: The basic functions gi{x) are calculated for x e [10^^,0.99] using numerical repre- 
sentations for the corresponding functions including the NIELSEN integrals. 

IRUN = 5 : The basic functions gi{x) are calculated for x e [10"'', 0.99] using the Mellin moment 
inversion by a complex contour integral. 

6.5 Subsidiary Routines 

A series of subsidiary routines representing mathematical functions and procedures are contained 
in the code. 



provides the value of the functions RES = (ZZ). Both ZZ and RES are COMPLEX* 16. 
6.5.2 /5«(z) 

SUBROUTINE BETA; (ZZ, RES) , k = 0..3 
provides the value of the functions RES = [5^^\tL), cf. Eq. Both ZZ and RES are COMPLEX* 16. 



6.5.1 V^('')(z) 



SUBROUTINE PSIfc(ZZ,RES) , k = 0..3 
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6.5.3 r(z) 

SUBROUTINE GAMMA (Z, RES) 

is the EULER Gamma function RES = r(Z) for complex argument, with 

r(z + i) ^ zr(z) . (141) 

Z and RES are C0MPLEX*16. For \z\ > 10, |arg|(2;) < tt the function is represented by 

log[r(^)] ~ z[\og{z)-l] + ^\og{27r/z) (142) 
1111 11 11 691 1 11 3617 1 

~l~ -» r» /-> /"V '} ~l~ -J r^ /-> /-A C t nntr\ <7 n /"v n /-> /"v 11 ~l~ 



12^ 360^3 1260^5 1680^7 360360^11 156 122400 ■ 
The SUBROUTINE GAMMAL(Z,RES) dehvers RES = log[r(Z)]. 

6.5.4 B(a,b) 

SUBROUTINE BETA (A, B, RES) 

is the EULER Beta function RES = i?(A,B) for complex argument, with 

A,B and RES are COMPLEX* 16. This function may be used to represent the Mellin transforms 
of the partonic input densities for complex argument. 

6.5.5 Harmonic Sums 

The code contains routines for the (alternating) harmonic sums up to depth 5 as REAL*8 
FUNCTIONS: 

DOUBLE PRECISION FUNCTION SUM1(I1,N) 

DOUBLE PRECISION FUNCTION SUM5 (II , 12 , 13 , 14, 15 ,N) 

The input parameters li may be either positive (non-alternating summation) or negative (al- 
ternating summation) integers. 

Example : 5'i_i,3(A^) is calculated by SUMS (1,-1, 3, N). 

6.5.6 Poly logarithms 
The REAL*8 FUNCTIONS 

DOUBLE PRECISION FUNCTION FLIi(x), i=2,3,4 
calculate the polylogarithms Lifc(x), /c = 2, 3, 4 for x e[— 1, 1]. 
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6.5.7 Nielsen Integral Si2(x) 

The Nielsen integral Si^2{x), cf. Eq. is calculated by 

DOUBLE PRECISION FUNCTION S12(x). 

6.5.8 Ii(x) 

The function Eq. ( [T3| ) is provided by 

DOUBLE PRECISION FUNCTION YIl(x). 

6.5.9 Numerical Integration 

REAL*8 FUNCTION DAIND(A,B,FUN,EPS,KEY,MAX,KOU,EST) 

The function DAIND yields the integral over the function FUN(x) from A to B. The integrand has 
to be declared as EXTERNAL in the calling routine. The corresponding algorithm was published in 
Ref. p]. We recommend to use the integrator setting KEY = 2. Here EPS denotes the demanded 
relative accuracy of the integral and MAX < 10000 the number of points at which FUN(x) is 
calculated by this adaptive integration. The output parameters KOU and EST refer to the number 
of points being used and the estimated accuracy reached. 

7 Summary 

We have calculated semi-analytic representations for the analytic continuations of the Mellin 
transforms of the set of basic functions through which the WiLSON coefficients and splitting 
functions which occur in hard scattering processes in massless field theories as QED and QCD 
can be expressed up to two-loop order. Here we aimed on high-precision representations which 
were performed using widely the algebraic and analytic relations of the Nielsen integrals and 
related functions being considered as well as their Mellin transforms limiting the necessary 
approximative representations to as few cases as possible. The expressions obtained are compared 
both to the Mellin moments at positive integer Mellin index with the values obtained by 
numerical integration and the representations in terms of harmonic sums. The Mellin inversion 
to X space provides a further test on the numerical accuracy of the expressions derived for the 
analytic continuations of the Mellin transforms for complex values. The FORTRAN-code ANCONT 
is provided. With the help of these representations the Mellin transforms for all two-loop 
quantities for the different space- and time-like hard scattering processes in the massless limit 
can be assembled. 

Acknowledgement. This work was supported in part by EU contract FMRX-CT98-0194 (DG 12 
- MIHT). 
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